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Abstract: The building block of the surface detector of the Pierre Auger Observatory is 
a water Cherenkov tank. The response to shower particles is simulated using a dedicated 
program based on GEANT4. To check the simulation chain, we compare the simulated signals 
produced by cosmic muons at various zenith angles with experimental data from a special 
Cherenkov detector equipped with a muon hodoscope. The signals from muon-decay electrons 
and the evolution of the charge with water level are also studied. 
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Introduction 

The surface array of the Pierre Auger Obser- 
vatory is a regular grid of cylindrical water 
Cherenkov tanks with a 1500 m spacing that 
samples the shower particles at the ground 
level [1]. The particles reaching the ground 
are mainly photons, electrons and muons, with 
mean energies around 10 MeV for photons 
and electrons and about 1 GeV for muons. 
Cherenkov radiation is emitted in the water 
by the electrons and the muons as well as by 
electrons produced by photons converted via 
Compton scattering and pair production. In 
each detector the water is highly purified and 
contained in a Tyvek® bag with a diffusively 
reflective white surface to maximize the path 
length of Cherenkov photons and thus their 
chance to be collected. The light is detected 
by three 9 inch photomultiplier tubes (PMTs) 
viewing the tank from the top. The PMT sig- 
nals are processed and digitised by 40 MHz 
Flash Analog-to-Digital Converters [1]. 
The full simulation chain of the tank response, 
from the produced light to the digitized signals, 
is performed with a dedicated program based 
on GEANT4 [2]. Here we describe the imple- 
mentation of the code and we will validate this 
tool by comparing the simulated response of 



an Auger tank to experimental data on muons 
crossing at different incident angles and water 
levels, and to electrons from muons stopping in 
the tank. 

Tank simulation framework 

The Auger tank simulation, which is a part of 
the Auger DPA Offline package [3], is based 
on the well established GEANT4 package. A 
dedicated module, called G4Fast, has been im- 
plemented to reduce the computing time. This 
module produces Cherenkov photons along the 
path of the injected particle and tracks them 
through the water until they are absorbed or 
they reach the active photocathode area of a 
PMT. The output is the number of photo- 
electrons as a function of time which is then 
processed by a different module simulating the 
PMTs and electronics response. 
Properties of a typical Auger water detector, 
such as the geometry of the tank and of PMTs, 
materials properties, etc.[l] are included in 
the simulation, but it is impractical to estab- 
lish with precision the details of the individ- 
ual properties of all the PMTs for each tank. 
Instead we use realistic average values given 
by their manufacturer: Photocathode Area = 
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426 cm 2 , Maximum Quantum Efficiency (QE) 
= 0.24, Collection Efficiency (CE) = 0.7. It 
must be noted that while these parameters di- 
rectly influence the number of photoelectrons 
produced by each PMT, the values are not cru- 
cial for understanding the responses of the de- 
tectors. Instead, the detectors are continuously 
calibrated with atmospheric muons: the mea- 
sured signals from showers are given in units 
of the charge of a vertical muon crossing the 
center of the tank [4] (Qvem or VEM). The 
simulated signals consequently also are given 
in VEM units. 

We also use the maximum water absorption 
length (L) = 100 m and the maximum Tyvek® 
reflectivity (R) = 0.940; these choices are dis- 
cussed below. 

Response of the tank for different 
water levels 

The water and Tyvek® parameters, L and 
i?, influence the propagation of the light in 
tank and hence how the signal decreases with 
time almost exponentially after the first reflec- 
tions. The chosen values of the parameters are 
such that the measured decay time is repro- 
duced well by simulations of simulated vertical 
muons. In fact, different pairs of values for L 
and R could reproduce equally well the exper- 
imental decay time. We present in this section 
a study meant not only to validate the simula- 
tion but also to disentangle L and R. 
An experiment was performed at the Auger 
site in a tank instrumented with scintillators 
to select vertical muons, where the water level 
was decreased over a week from 120 cm to 
75 cm. The charge deposited by the vertical 
muons was measured for every 2 cm drop in 
water level. Simulations were performed using 
G4Fast with the water/Tyvek® parameters L, 
R given in the previous section and with an- 
other set, where the Tyvek® quality was im- 
proved and the water made more attenuating 
(L = 30 m and R = 0.973). 
Working with different water levels changes the 
relative influence of water and Tyvek®. Less 
water means a reduced volume where photons 
impact more often on the Tyvek®: the impor- 
tance of the reflectivity is expected to increase 
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Figure 1: VEM charge relative to the track- 
length as a function of water level 

as the water level decreases. The main effect 
of less water is to decrease the track length 
of vertical muons, and so the deposited charge. 
However, at the same time, the charge per unit 
of length increases, and we expect an enhanced 
effect for the larger R. The VEM charge, nor- 
malized to the tracklength, versus water level 
is shown in Fig 1 for data and for simulation 
with the two different sets of parameters. The 
expected effect of Tyvek® is clearly demon- 
strated and the data are completely consistent 
with the chosen parameters. 

Response of the tank to vertical 
and omni-directional muons 

The basic calibration information is the charge 
deposited by vertical and central throughgoing 
muons. In this section, we compare the simu- 
lation with vertical muon data as a first test of 
the simulation. 

The water tank, in its normal configuration, 
has no way to select only vertical muons: how- 
ever the distribution of charges deposited by 
omni-directional muons has a peak which is 
well correlated with the VEM charge [4]. The 
peak is at 1.09 VEM, measured in an Auger 
tank instrumented with a muon hodoscope [5] . 
This ratio is an essential parameter to be repro- 
duced by the simulations. Using G4Fast, we 
simulate vertical muons (i.e. passing through 
the center of the tank and crossing the entire 
volume of water) as well as omni-directional 
muons. We use a realistic spectrum of multi- 
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directional electrons, muons and photons with 
the energy spectrum from reference [6]. We 
assume, as an approximation, the same zenith 
angle distribution, f(6) — cos 2 (#) sin(#), for all 
particle types. We use a low-threshold trigger 
requiring a 3-fold coincidence over 0.15 iy^ M 
in each PMT as in real data (being I^^ k M the 
average of the peak in the pulse produced by 
vertical muons). Fig. 2 shows the compari- 
son of the experimental and simulated charge 
distribution for the omni-directional muons in 
units of VEM: the simulation reproduces the 
data well. The position of the peak is found to 
be 1.09 VEM, as in data. 




Signal(VEM) 



Figure 2: Charge distribution (VEM units) 
for omni-directional particles. Full line corre- 
sponds to data and dashed one to simulation. 

Response of the tank to inclined 
muons 

To validate the simulations of inclined muons 
we compared the simulated values of the 
recorded charge with measurements from a 
test tank in Orsay, similar to the Auger ones, 
where signals from atmospheric muons could 
be recorded at different zenith angles. Two 
movable scintillators were placed at the side 
of the tank, triggering on muons arriving with 
different angles, disentangled from correlated 
shower events by means of proper timing [5] . 
The scintillators were located for each incident 
angle in two opposite positions as shown in 
Fig. 3. As the average energy of the muons 
increases with zenith angle, we use the energy 
parameterisation given in [6]. The simulated 
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Figure 3: Set-up of the Orsay tank. 

and experimental charges are plotted versus 
the muon tracklength in Fig. 4. A deviation 
from linear behaviour is observed in the data 
as the zenith angle increases, due to the ap- 
pearance of direct non-reflected light, in par- 
ticular when muons cross close to the PMT. 
This behavior is reproduced well (within 10%) 
by the simulation. 
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Figure 4: Mean charge of the 3 PMTs as a 
function of tracklength in tank (both quantities 
are normalized to the VEM) 

Response of the tank to electrons 
from muon decay 

Muons decaying in the tank produce electrons 
with a well known energy spectrum, the Michel 
spectrum, with an end point at 53 MeV and an 
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average value of 37 MeV. The measurement of 
the Cherenkov light produced by Michel elec- 
trons provides a reference point for the tank 
response to low energy electrons. A dedicated 
trigger was implemented to select muon decay 
events and allowed measurement of the ratio 
between the Michel spectrum peak and that 
of vertical throughgoing muons: the average 
value for the 230 tanks in operation at that 
time was found to be 0.13 with a tank-to-tank 
spread of 0.01. A study with G4Fast used the 
same algorithm as used in the data to select 
the decaying muons. Crossing muons were gen- 
erated with an appropriate angular distribu- 
tion and energy spectrum [6]. The simulated 
ratio between the electron and muon charge, 
0.13 ±0.01, agrees well with the measured one. 



-1.3 



£ 1.2 



1.1 



0.9 



• G4 SIMULATION 



. t DATA 



5.5 6 6.5 7 7.5 8 

mu/e ratio 

Figure 5: Evolution of muon to electron charge 
with the water level 

Michel electrons are absorbed in less than 25 
cm and therefore are almost insensitive to a 
change in the water level, while vertical muons 
produce a signal mostly proportional to it. A 
linear dependence of the ratio between elec- 
trons and muons versus water level is expected. 
An experiment was performed to test the influ- 
ence of water loss: a tank was slowly drained 
starting from its normal level at 1.2 m and fin- 
ishing at 0.895 m [7]. Four ratios were ob- 
tained during this process, shown in Fig. 5 as 
triangles. Simulations were carried out for the 
same water levels and the e/fj, ratios computed 
as described previously. The compatibility be- 
tween the simulation points and the data is 
visible in the same figure. 



Conclusions 

The simulation of the Auger water Cherenkov 
tank is accomplished by a module based on 
GEANT4, designed to reduce the computing 
time (so called G4Fast). We have described a 
variety of tests of the simulation versus data: 

- Vertical muons. The ratio between the VEM 
charge (the basis of the calibration of the 
Auger surface detector) and the average charge 
detected for omnidirectional muons has been 
measured in a dedicated experiment to be 0.92. 
The same value is found using G4Fast sim- 
ulations of both vertical and multidirectional 
muons, and their spectrum also is reproduced 
well. 

- Inclined muons. Due to their increased path- 
lengths in the tank, inclined muons yield a 
larger charge. The behaviour of the signal ver- 
sus tracklength has been measured in an ad hoc 
experiment: the simulated charge response for 
different muon directions is well represented by 
G4Fast, including the effects of direct light on 
the PMTs. 

- Electrons from muon decay. We have mea- 
sured and simulated the ratio between the 
charge peak from Michel electrons and the 
VEM peak. The average experimental ratio 
(0.13, with a tank-to-tank dispersion of 0.01) 
is reproduced well by G4Fast. Good agreement 
is found as well for the evolution of this ratio 
with changes in water level. 
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